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Abstract 

It is possible to determine an off-shell propagator for heavy quarks to order 
1/m in mass and in any covariant gauge £, which applies universally to all the 
quarks, by using the gauge technique. The result for the leading behaviour of 
the propagator is 

where v is the quark velocity, A is a QCD mass scale and = a(2 + ^)/37r. It 

is totally reliable in the infrared limit and accounts for soft-gluon corrections 

to the fermion in internal loops. 
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*in memory of Abdus Salam, to whom we owe the gauge technique. 
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I. INTRODUCTION 



The QCD quark Lagrangian (N flavors) is endowed with a higher symmetry in the limit 
of equal quark velocity which applies even when the quark masses are different JIJ. Thus 
it generalizes the old U(N) x U(N) supersymmetry |2| for the equal mass case. Provided 
that the momentum transfer to the gluons is not much greater than the QCD scale one can 
thereby deduce a number of relations between transition amplitudes, which seem to be borne 
out by experiment. It has become customary || to attribute a velocity v to the constituent 
heavy quark so that the momentum of the quark in a bound state is written p = mv + k, 
where k denotes the residual quark momentum, itself associated with the light material that 
makes up the hadron. As a result one can show that the "free" quark propagator in the 
m — > oo limit is simply given by S = (1 + 7 ■ v)/2v ■ k and one can use this in subsequent 
leading order calculations of various matrix elements. 

In this paper we would like to show that one can improve on S by taking account of soft 
gluon corrections. The result of the dressing is to provide a propagator which contains the 
characteristic QCD scale A and which coincides with the free one in the limit of vanishing 
gluon coupling a = g 2 /^. Thus this propagator applies just as well to all the quarks in the 
heavy mass limit [[| and does not jeopardise the prevailing higher symmetry. In order to 
derive it we use the gauge technique || for QCD, which is known to be a reliable method in 
the infrared and ultraviolet limit ||. The technique produces a self-consistent equation for 
the quark spectral function in any gauge, from which the propagator follows [0]. In the next 
section we set out the velocity projector decomposition of the propagator. Next we derive 
the effective vertex for soft gluons and finally we solve the equation in question, obtaining 
the result quoted in the abstract; there we also compare the result with QED where the 
scale A is missing. 
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II. VELOCITY PROJECTIONS 



When one substitutes p = mv + k in the free quark propagator, the resulting expression, 

1 



S(P) 



(1) 



777.(7 • v — 1) + j • k' 

has to be taken in the limit m — > 00 in order to discern the resulting (leading) dependence 
on four- velocity v (with v 2 = 1.) Now any quark matrix M for a particular flavor can be 
decomposed into projections using P± = (1 ± 7 • v)/2 according to 

M = P. M. . I'. + I'.M. P + P M . I'. + P M P . 

This has the effect of resolving the 4x4 matrix into four separate 2x2 matrices: 



M 



M_ . M_ 



In this basis, 



V " 1 / 



7/, 



In particular the inverse free propagator, resolves to 

S^ip) = m (7 • f — 1) + 7 • 
and correspondingly, 

Sip) 



( \ 

v„ 7m 

v 7m 



ft., 7 .t ^ 

7 • A; —A; • t> — 2m 



^ k ■ v + 2m 7 • A; 



k 2 + (k ■ v) 2 + 2mA; • i> 



V 



7 • A; —A; • t> 



Thus, up to order 1/m 2 , the free S decomposes into 



Sip) 



I y k 2 7-fc \ 



A; • v 



2mk-v 2m 



7-fc 
\ 2m 



k-v 
'2m J 



m z 



from which one infers that the leading large component is S ++ ~ 1/k ■ v. 



(2) 



(3) 



(4) 



Let us now consider the full propagator in a covariant gauge, which is best written in 
the Lehmann-Kallen spectral form for our purposes: 

p(W) dW 



S(p) 



dW = 



—m poo 
+ 

oo 



)dW } 



(5) 



7 • p - W + ir}e(W) ' 

for an ordinary sort of particle. In the heavy quark limit we anticipate that the negative 
energy cut is 'far away' and that the main contribution from soft gluons will arise in the 
vicinity of W = m, to within a region of order A. The free propagator is of course obtained 
just by setting p(W) = 5(W — m) above. Put W = m + lj, p = mv + k and take velocity 
projections as in (4), to obtain 



S(p) 



p(cj) d/jj 



2m(v ■ k — lo) + k 2 — oo 2 



2m + k ■ v + lj 7 ■ k 
7 ■ k —k ■ v + u) 



(6) 



We see that the dressed propagator is still dominated by its S ++ component, which assumes 
a very simple form, / dujp{uj) /{v ■ k — lj) despite the inclusion of QCD interactions. More 
generally, to order 1/m, we get 

/ - . (vkf-k 2 T k \ 



S 



p(cj) du 
v ■ k — LJ 



1 + 



7-fc 



V 



2in(v-k—ui) 2m 
LO—v-k 



7-fc 
2^7 



2m / 



III. APPLICATION OF THE GAUGE TECHNIQUE 



The next stage involves solution of the Dyson- Schwinger (D-S) equation for the propa- 
gator, while taking cognizance of the (longitudinal) Ward-Takahashi identity, 



(p - p , rS(p')T^p',p)S(p) = S(p') - Sip). 



(7) 



because it can lead to a self-consistent equation for S, if we ignore certain transverse terms 
in the vertex T M . If one were to incorporate the transverse Takahashi identity 0] as well, 
one would in effect be solving the full field theory; but the transverse identity — in anything 
but two dimensions || — brings in other vertices leading to a system of equations which is 
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actually not closed, unless one makes further drastic truncations [10]. Alternatively if one 



knew the full solution of the D-S equation in any particular gauge, one would be able to 



determine it in any other gauge via the Landau-Khalatnikov-Zumino gauge covariance |TT 
relations. 

We have none of these luxuries. The gauge technique does its best to solve the equation 
(7) in the form stated, while making sure that the singularities in the non-truncated Green 
function are properly included. It does not solve the inverse form of the equation (7) because 
that would give a linear relation between T and S and produce a difficult nonlinear 
equation for the inverse propagator; besides which, it is not obvious how to handle the 
heavy quark limit for the inverse propagator — which is dominated by its Szh projection, 
conversely to S ++ . The gauge technique starts off with the obvious solution to the vertex, 

S(p')Tl(p',p)S(p) = f dWpiW) i-— 7, (8) 

as a weighted mass integral. One can readily check that (8) obeys (7) automatically, but 
of course (8) is subject to transverse additions r ; these are unknown unless one has some 
knowledge about them through perturbation theory |12 or examines equations for higher- 



order Green functions |L3| or makes use of the transverse identity ||, which is essentially 
equivalent. It is worth pointing out that such transverse terms are soft, vanishing with the 
vector meson momentum. For that reason the gauge technique is a clearly reliable tool in the 
infrared limit, though it is also gauge- covariant in the ultraviolet regime || as it happens; 
it is only at intermediate energies that transverse corrections to (8) play an important role. 

Returning to heavy quarks, let us expand the solution (8) in powers of 1/m by writing 
it as 

S(p')Tl(p',p)S(p) = [ dup(u)— -i- 1 ,— -i- (9) 

p J m(7 ■ v — 1) + 7 ■ k — uj m(7 • v — 1) + 7 • k — u 

and taking velocity projections: 



p(uj) dui 



2m(v-k'-uj) 2m(v-k-oj) 



(v ■ k' — 10) (v ■ k — 00) 



^ 2m 2m 



y-k' ■v^+y^m-k' -v) 
2m 



(10) 



up to order 1/m 2 . Note that representations (8) and (10) are exact for small p — p' in the 
complete Green function (7), for the reasons we have already given. 

The next step is to use the approximation (8) in the D-S equation, which we write in 
the renormalized form, 



Z-i = ( 7 • p - m + Sm)S(p) + / d 4 qS(p)T,(p,p - q)S{p) lv D^{q)-. (11) 

Recalling the connection, Z~ x = J p(W)dW, the spectral form of the equation is, 

p(W) dW. 



7 • p — W 



[W - m + 8m + E(p, W)} = 0, 



(12) 



where 



(2tt)< 



(13) 



is the self-energy for a quark of mass W due to gluons in first order perturbation theory. 

At this point we carry out the heavy quark expansion and take velocity projections to 
arrive at 



/ 



p{u) dw 
v ■ k — uj 



( 



(u + 8m) 



1 + 



2m(v-k—uj) 2m 



■y-k 
2m 



LO — V-k 

2m J 



£-+(i> • k, uj) E (v-k,uj) J 



= 0, 
(14) 



where, after summing over colours (hence the factor of 4/3),, 

/ 



E(t> • k, uj) 



4g 2 /3 f d 4 qD^(q) 



(2tt) 4 J v-{k-q)-uj 
and up to order (1/m). 



v»v v [1 + 0(1/ m)] v^ v + 0(l/m) 

\ 2m 2m 2m 2m 



(15) 



IV. THE SPECTRAL EQUATION 



To make any further progress and determine the spectral function p and thence the 
propagator, we need to make some further approximations/assumptions about the behaviour 
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of the gluon. It is generally accepted that the gluons are massless so that the propagator 
D(q) is at least as singular as 1/q 2 ; it is also known that in the ultraviolet regime this is 
subject to well-defined logarthmic damping; the behaviour for small q 2 , where the strong 
force enslaves colour, is more mysterious and there have been suggestions that D(q) could 
be as singular as 1/g 4 , that it plateaus or even that one should not be using QCD at all but 
an effective field theory incorporating chiral symmetry with real mesons. What is certain is 
the occurrence of a mass scale A demarcating the ultraviolet from the infrared regime of D. 
As we are only interested in soft gluon effects on the heavy quark lines, we will adopt a gluon 
propagator which implies masslessness, which cuts off in the ultraviolet and which introduces 
the fundamental QCD mass scale. For our purposes it is enough to use an effective 

D^(q) = (-V„u + i^) q2 ^_ q2 y (16) 

knowing its limitations full well. It incorporates the main things we want and also includes 
a covariant gauge parameter £. If other readers wish to modify D with a more sophisticated 
and perhaps more realistic expression, they can repeat our calculations below; while that is 
sure to alter the precise form of our answers, we believe it will not affect the main features 
of our results in a very significant way. 

Returning to the largest component of (15), we have to consider the spectral equation, 

0= f P ^ du [u + 8m + i: ++ (vk,u)], (17) 
J v ■ k — uj 

■ k ^ )=l j2^i d q v.(k-q)-U ~ q^ 2 ) ■ (18) 

A straightforward but messy calculation gives 



E ++ (u;,u>) = — 



{u - u')(i + 1 - r A2 ; ) in v A ' + -a^ - ^')+ 



1 (u - UJ') 2 I LU-LU' - J(UJ - UO') 2 - A 2 

=- 1 - e A , ] ) v /(cu-^) 2 -A 2 ln V , ^ 



(19) 



to leading order and in a general gauge. First we solve the spectral equation (17) 
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J ad p{uj') 



(20) 



in the Fermi-Feynman gauge £ = to discover what is going on. By taking the imaginary 
part of (20), we obtain 



oo 

/\ 7 / 



= (uu + 5m - 2ak/ 3) p(u) + (4a /3ti) / p(d) du . (21) 
This has the solution 

oc (w + 5m - 2aA/3)' 1 ~ 4Q/3w 

and, since the self-mass is given by 

J (u + Sm)p(u) dcu = 0, 

this fixes 5m = 2a A/ 3. It makes good sense, being governed by the QCD mass scale and 
gluon coupling. One last matter is the proportionality factor: we must ensure that p(ui) 
reduces to 5(u) when a — >• 0. Hence we choose the overall constant so that the result for 
the spectral function for £ = is neat and compact, namely 

Pt=o(u>)= >TV t~~7^~t ( ~ ) • (22) 



ujT(-4a/37r) \co 

A bonus of this choice is that the heavy quark propagator simplifies to the elegant non- 
perturbative expression, 

s ^.. ) = r ( i + -)^-I r (-— ) . (23) 

In the limit a — > one recovers the free result l/(v ■ k) for S. 

Now we turn to the general gauge £. Noting that the ^-dependent part of S++(f • 
vanishes at the threshold v ■ k — cu, the gauge- dependence of p arises purely from the 
imaginary part of In this way (21) gets modified to 



4a r°° 

= (uj + 5m - 2a/3n)p(uj) + — / dd p{uj') 

07T Jul 



i >-u/) 2 



1 + 



2 * A 2 



(24) 



Once again the self-mass condition requires 6m = 2a A/3, which is satisfyingly gauge- 
independent. The resulting integral equation for the spectral function is a little bit harder 
to solve now; nevertheless one may establish that it reduces to a generalised hypergeometric 
function: 

1 /A\- 2a « „ A _ ,_ a£u 2 \ 



" (w) = ^rF^y .^i + °e.s/2 + <*;-jj^J, (25) 



where ctg = a (2 + £)/37r. Thereupon [14] the heavy quark propagator becomes 



% • k) = r(l + 2a,)i±^ (-^) F 2 (l + Of, 3/2 + a,; + 0(l/m). 

(26) 

This is universal to all the quarks and could be used to estimate the gluon corrections in 
loops which result from dressing fermion lines and their vertices. However, a word of caution: 
the result (26) does not take account of gluon self-interactions; those will somehow need to 
be included separately in heavy quark calculations. 

There is one further test of our work. One needs to verify that to order 1/m the other, 
'small component' sectors in the propagator velocity projections are correctly determined 
by (25), because they are fixed in terms of the leading p. We have indeed checked this out: 
the S.\ sector produces precisely the same equation as (24), while the S sector is nothing 



but the self-mass condition, J(u + 5m)p{ui) du = 0, which we have already settled [|T^ . 

Finally it is worth comparing the answers with the gauge technique solutions for scalar 
QED say. Those solutions do not have the benefit of an intrinsic cut-off; rather the source 
mass itself acts as the cutoff and the results read, 

(W 2 lm 2 - l^r 1 - 2 ^ w 2 

m 2 p{W) = ± ' \ 2 Fi{-a^ 1 - a ? ; -2a 5 ; 1 - — ); a ? = (f + 2)a/47r 

T(— 2a^) m 2 

m 2 S(p) = r(l + a e )T(2 + 0^2^(1 + ^,2 + ^; 2; p 2 /m 2 ) 

We notice a strong similarity with (25) and (26), which becomes greater when one substitutes 
p = mv and expands to order 1/m. The only change is that m takes the place of A as 



the argument of the hypergeometric function. It only remains to obtain the non-leading 
behaviour of S in the various sectors. This entails solving the spectral equation up to order 
1/m and is where our use of the approximations (9) and (16) start to look a bit suspect 
because they are connected with gluons which carry off appreciable momentum. It is a 
nice subject for future research since it portrays the mass-dependence of the heavy quark 
Lagrangian. 
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